Strongly correlated hopping and many-body bound states 
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We study a system in which the quantum dynamics of electrons depend on the particle density in 
their neighborhood. For any on-site repulsive interaction, we show that the exact two-body and 
three-body ground states are bound states. We also discuss the finite density case in a mean- 
field framework and we show that the system can undergo an unusual transition from an effective 
attractive interaction to a repulsive one, when varying the electron density. 
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I. INTRODUCTION 

Correlated hopping models have been the subject of 
many studies within various contexts. First proposed by 
M. E. Foglio and L. M. Falicov to decribe mixed valence 
solids jl], they have also been widely used to study the 
organic conductors |p|-p[ in order to take into account 
bond-charge effect ff. Finally, these models where the 
probability of an electron to move depends on the par- 
ticle density, have been proposed to mimick effective at- 
tractive interaction between electrons high- Tc supercon- 
ductivity 1^.^) and have provided rich phase diagrams 

In a completely different framework, we have recently 
described a localization phenomenon induced by the 
magnetic field that occurs for special geometries and for 
special values of the magnetic flux |27| . This surprising 
effect has been experimentally observed in superconduct- 
ing wire networks [^8| , and in two-dimensional electron 
gas . We have also studied the influence of electron- 
electron interactions on such systems and we have shown 
that a Hubbard-like term (on-site repulsion) was able to 
delocalize two particles initially confined in a given so- 
called Aharonov-Bohm cage [ ^J3l| . These results have 
led us to formulate the simple toy model presented here 
in which we introduce this delocalization process directly 
in the hamiltonian by imposing that an electron can move 
only if another electron is in its close neighbourhood. Of 
course, we do not claim to capture all the physics of the 
interacting Aharonov-Bohm cages with this simple one- 
dimensional system, but we think that it can help us in 
understanding the delocalization process induced by (re- 
pulsive) interactions. 

This paper is organized in two main parts. In the first 
one, we exactly solve the two-body (section^ and three- 
body (section HI) problems and we show that the ground 
state is always a bound state for any strength of the on- 



site repulsion. Therefore, we give a simple picture in 
terms of a graph in the Hilbert space that makes the 
physical interpretation clearer. In the second part (sec- 
tion IV), we propose a mean-field like approach for the 



finite density case by considering the low-energy excita- 
tions above the Fermi sea. We show that in the vicinity 
of a curve int its two parameter space, the initial hamil- 
tonian can be mapped onto an effective Hubbard model 
with an interacting term that depends on the particle 
density and that can be either repulsive or attractive. 

We consider interacting spin 1/2 fermions in a one- 
dimensional chain of linear size L = Na with peri- 
odic boundary conditions described by the following 
hamiltonian : 

1,(7 

+[/^ni,|n,,x, (1) 



where (resp. Ci^a) denotes the creation (resp. annihi- 
lation) operator of a fermion with spin a, rii^^ ^ cj (j Ci,o- 
the density of spin a =t, i fermion on site i, and (...) 
stands for nearest neighbor pairs. The kinetic part of the 
hamiltonian (^ allows a particle of spin a located on a 
site i to jump on a neighbouring site j only if there is 
already a particle either on site i or on site j. For simpli- 
city, we restrict our analysis to the repulsive case U > 0. 
However, since the structure is bipartite, the spectrum 
of H is odd under the tranformation U ^ —U iBSl. 



II. THE TWO-BODY PROBLEM 

The single-body problem is trivially solvable for this 
model since the particle can neither move, nor interact. 
In this case, the spectrum consists in one eigenvalue e = 



1 



which is A''-fold degenerate. So, let us pay attention to 
the two-body problem. Denoting \i) the orbital localized 
on the site i, the two-body state space is generated by : 



(2) 



We only consider here the orbital degrees of freedom 
since for two particles, the spin degrees of freedom are 
completely determined by the symmetry of the orbital 
wavefunction. This problem can be mapped on that of 
a single particle moving in the graph displayed in Fig. |l| 
where each site of the lattice (with appropriate 
boundary conditions) corresponds to the ket \i,j). Us- 
ing the invariance of the system under a translation of 
the center of mass (direction (1,1) in Z^), diagonalizing 
H can be simply achieved by introducing the following 
Bloch waves : 
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(3) 
(4) 



where b = a%/2 and where K = 2ti j/Nb [j e [0, iV - 1]) 
is the total momentum of the two particles. 
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FIG. 1. Representation of the two-body problem as a graph 
standing in the lattice and its reduction after projection 
along (1,1). tK =t{l + e-''^''). 



Indeed, one has : 
{ipo{K)\H\^i{K')) 



-til 
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(5) 



(6) 



where 5k. k' is the usual Kronecker symbol, so that the 
eigenvalues of H are given by £ = for the triplet states 
(non sensitive to U) and for all the trivial configurations 
corresponding to isolated particles, and by : 



e±{K) ^\{y± y^U^ + 32t^ coa^ {Kb /2) 



(7) 



for the singlet states. Note that this result has already 
been obtained by Hirsch in a different context p^ , ^ . 
The most surprising fact is the emergence of a dispersive 
band S-{K) < associated to two-body bound states 
for any U. In particular, the ground state is obtained for 



e-{K = 0). This can be understood invoking the com- 
petition between the kinetic term that lowers the energy 
and that is enhanced when particles are close together, 
and the interaction term that favors the opposite situa- 
tion. An interesting issue is to know whether this feature 
still holds for more than two particles. To tackle this 
task, we shall now analyze the three-body problem that 
is, fortunately, still tractable analytically. 



III. THE THREE-BODY PROBLEM 

As previously, we consider the three-body problem in 
the one-dimensional chain as a single-body problem in a 
graph standing in the lattice. Using, once again, the 
invariance of the system under a translation of the center 
of mass (direction (1,1,1) in Z'^), it is convenient to build 
the following Bloch waves : 
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+ 71, j + n, fc + n), (8) 



where c = a-\/3 and where K — 2-1:] /Nc (j € [0, iV — 1]) 
is the total momentum of the three fermions. Remark 
that \(pij,k) oc \ipi'j'M') if ihj.k) = {i' , j' ,k') mod {1,1, 1) 
so that one must only consider non equivalent Bloch func- 
tion. Here, we have chosen to affect the phase 1 to the 
ket \i,j,k) such that (i^+j^ + fc^) is minimum. After pro- 
jection along (1,1,1) the graph is the so-called sextopus 
displayed in Fig. |[ 




FIG. 2. Representation of the three-body problem graph 
after projection perpendicularly to the direction (1,1,1) in 
the lattice. The sites denoted ijk correspond to the ket 
\'P^.J.k{K)) und z = -2te-"^\ 

Note that in this latter graph, we have not put the site 
corresponding to \(poflfi{K)) since it is forbidden for spin 
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1/2 fermions but that would be allowed for bosons. It is 
readily seen with this representation that there are two 
types of regions : (i) the "bulk" where the three particles 
are close together and where there could possibly exist 
three-body bound states ; (ii) the "legs" corresponding 
to a situation where one particle is motionless and the 
two others propagate together. One also observes that 
K plays the role of a "pseudo-magnetic field" for the 
system so that the ground states will be obtained for 
K = (zero field condition). For an arbitrary value of 
K ^ mod tt/c, the sextopus has the dihedral symmetry 
Z?3 and can be easily diagonalized for any N. Note that 
when the three electrons have the same polarization, the 
system is frozen since the particles cannot move. It im- 
plies that the energy of the quadruplet 5" = 3/2 is simply 
0. The only interesting situation thus arises in the sector 
S' = 1/2 whose spectrum is shown in Fig. || as a function 
of K for a given U. 
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FIG. 3. Spectrum of the three-body problem as a function 
of Kc for (7 = 1 and TV = 100 {t = 1) . 

One clearly observes two distinct components. First, 
for a fixed K, there are two dispersive bands that corre- 
ponds to scattering states (one two-body bound state -t- 
one motionless particle) propagating ballistically in the 
"legs". Their precise form is given by the equation (|^). 
Second, one observes much more interesting states out of 
those bands. As we shall see thereafter, these states are 
actually three-body bound states and the most remark- 
able fact is that the ground state is always given by such 
a state for any U. To analyze more precisely this surpris- 
ing phenomenon, we focus on the case K = for which 
the ground state is obtained. In this case, the symmetry 
of the sextopus is Dq and it is possible to look for bound 
states in each representation schematized in Fig. |[ Each 
site represented in Fig. ^ symbolizes a state vector with 
fixed angular momentum. Note that the representations 
indexed by ^ = 0, 3 corresponds to a fully space symmet- 
ric wavefunction that can, in no case, be an eigenstate 
for three spin 1/2 fermions. 
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FIG. 4. Non trivial graph of each irreducible representation 
of De indexed by I G [0, 5]. Ui and ti are defined in the text. 

We thus look for a bound eigenstate j^*) = J2n''l^ri\n) 
where \n) is the ket corresponding to the n*'' site of 
the half-chains labelled by Z £ [1,2] (or equivalently 
I e [4,5]). For aU n > 3, we set Vn = CAe-^("-3)6 
for n odd and t/jn — D\e~^'^'^^^^^ for n even, and we seek 
for complex A such that 3?(A6) > and 3(A6) = modir, 
this latter condition ensuring to have a real eigenenergy. 
Such a state is an eigenstate if the secular equations are 
satisfied : 

(4(A)-C/,)Vi --2t, ^2 

ef(A)V2 = -2t,(V'i+V'3) (9) 
(e|(A) - U) V^3 = -2t, ^2 - V2t ^4, 
where Ui = U — 4t cos{2ttI / 6) and ti — i cos(7rZ/6) and : 

e^i^X) ^^(^U± ^Ju^ + 32t^ cosh^(A6/2)^ . (10) 

We emphasize that the periodic boundary condition adds 
restrictions on A but since the amplitude decreases ex- 
ponentially, they are not relevant provided N is larger 
than the localization length 1/A. A simple inspection of 
the secular system (^ allows to show that it always has, 
at least, one solution Aq (for any U > 0) associated to 
e^{Xo) which is, thereby, the ground state of the system. 
We have plotted in Fig. gthe dependence of the two-body 
and three-body ground states energy as a function of U 
for spin 1/2 fermions. 
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FIG. 5. Variation of the ground state energy of the 
two-body (□) and three-body (•) problems as a function of 
the interaction U {t = 1). 

When U goes to infinity, both energies converges to- 
ward zero but, as noticed previously, the three-body 
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ground state has always a lower energy than the two- 
body one. In addition, the localization length I/Aq tends 
to infinity which sounds quite natural for strong repul- 
sion. 

IV. THE FINITE DENSITY CASE 

A natural question then arises : are there still many- 
body bound states at higher densities ? Indeed, if we now 
consider four electrons, we can figure out more complex 
scattering processes where two two-body bound states 
collide with each other, or where a bound pair oscil- 
lates between two isolated particles. In these cases, one 
may expect the emergence of four-body bound states. 
Unfortunately, this problem is too complicated to be 
completely analyzed by elementary methods. We have 
searched for the possibility of binding together two two- 
body ground states by means of a variational method. 
In this approach, the four-particle state space has been 
decomposed in two sectors. The first one corresponds 
to two isolated two-particle states whose centers of mass 
are separated by at least 2.5 lattice spacings. In this 
region, we assumed an exponential decay of the wave- 
function with the relative separation of the two center of 
mass positions. The second region corresponds to real 
space configurations where each of the four particles is 
close to any other. In the simplest version, this defines, 
after taking global translation invariance into account, a 
finite set of 18 configurations. It is important to note 
that two two-body bound states collision may generate 
real space configurations where a bound pair escapes to 
infinity leaving two unpaired electrons which may remain 
at an arbitrary large final distance from each other. This 
is illustrated on the sequence shown in Fig. |^. 

t 1 t I 

t tl i 

t ■ i tt ■ ■ 

t ■ i t I ■ ■ 

t ti I 

t ■ I tt ■ 

t ■ ■ I t I ■ 

t I tl 

FIG. 6. A possible way to obtain an isolated pair and two 
frozen electrons from a four-body compact configuration. 

However, the typical decay length for the escape of a 
single pair is likely to be significantly shorter than the 
one for two bound pairs. Thus, we have not kept these 
asymptotic states with a single pair in our variational 



approach. This search for four-particle bound states has 
always failed so far. Of course, this lack of evidence does 
not constitute a proof of the absence of such bound states 
in the spectrum, but we believe that four-particle bind- 
ing is unlikely in this model. The main reason for that is 
the Pauli principle which severely restricts the possibil- 
ity of particle hopping in a close packing configuration. 
Therefore, we suggest that the finite but small density 
system will form a Luttinger liquid of bound pairs. The 
absence of four-particle binding may mean that residual 
interactions between these bound pairs are repulsive, al- 
lowing a well-defined thermodynamic limit. This picture 
of a fluid of bound pairs is in fact reminiscent of the at- 
tractive Hubbard model js^. Indeed, a detailed analysis 
of the Bethe ansatz spectrum shows that this model 
is precisely described in its low-energy limit in terms of a 
Luttinger liquid of spinless two-particle bound states. By 
analogy, we thus also expect a gap in the spin excitation 
spectrum of our model. Actually, the ground state and 
the thermodynamical properties of both models should 
be quite similar. However, dynamical quantities might 
produce some meaningful differences. For instance, it 
would be very interesting to investigate in more details 
the behaviour of the single electron Green's function. In 
our simple-minded picture, adding an electron to a state 
with an even number of particles indeed leaves an un- 
paired spin. By contrast to the usual spin charge sep- 
arated liquid and to the attractive Hubbard model in 
particular, our model opens the possibility of sustain- 
ing three-body bound states between this extra electron 
and a pair taken from the Luttinger liquid ground state. 
This interesting question clearly deserves further stud- 
ies. Instead of adressing this issue here, we shall turn 
to a very simple description of the finite density system. 
For large particle densities, the two-particle bound states 
have to overlap in real space, and it is no longer clear 
that they provide a good basis to understand the ground 
state properties. In the following, we shall assume that 
the quantum state of the finite density system is not too 
far from a Slater determinant of plane waves. Let us now 
write the hamiltonian in momentum space : 

k,k' .,p,p' .a 

f{p,p'; k\ k) cl^^cl, _^^Ck'.~aCk,a, (11) 

where Cp^a = J2n e™^°c„,CT and p is the momentum 
of the state. We keep the same notation c for operators 
either in real or in momentum space since no confusion 
is possible here. The interaction function is given by : 

/(p,p'; k', k) - t[cos{pa) + cos(p'a) + 

cos(fc'a) -I- cos(A:a)]. (12) 

The expectation value of H taken on any free particle 
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state written using the plane wave single particle basis is 
given by : 



(13) 



k,k' 



To determine the most stable Slater determinant of plane 
waves, we consider a single particle-hole excitation away 
from a reference state whose occupation numbers in mo- 
mentum space are denoted by nfl^. For instance, we 

consider : Uk,]- = n^°| -I- Sk,p - Sk,q and nk,i = The 
average energy change induced by this particle-hole ex- 
citation is : 



{SH)p^g = - ^[/(p, k'; k',p) - f{q, k'- fc', q)]nf^ (14) 



-4t— [cos(pa) — cos(ga)]. 



(15) 



where iVo- is the total number of spin a electrons. This 
shows that the average energy is minimal provided the 
occupied state (momenta) fills the usual Fermi intervals 
] — kp.a, kp^a] for a =t, i- From these considerations, and 
from the expression ([III), it is very natural to separate 
the diagonal and off-diagonal part of H . We thus have : 



U 
N 



k,(T 



teff,aCOs{ka)cl^^Ck 



(16) 



with teff,cr = 2tN-„/N. Up to a global shift, Hdiag ex- 
hibits the same structure as a pure hopping hamiltonian 
for free particles, but its hopping term depends on the 
particle density as it is induced by a two-particle inter- 
action process. Along the same line, one has : 

Hoff-diag = ^ 5k+k' .p+p' (1 — Sk,p) X 

k,k' ,p,p' ,(T 

f{P,p',k',p) cl^^cl, _^Ck'-aCk,a- (17) 

Let us assume that the ground state is not too remote 
from the usual non-interacting Fermi sea. In this case, 
the four wave vectors involved in (^) are close to Fermi 
points so that Hoff-diag is well approximated by : 



off-diag 



— ^ 6k+k',p+p' - 6k,p) 



k,k' ,p.p' ,(7 



u, 



5£ t t 

r) ^p.a'^p' ,-(y^k' ,~iyCk.o 



(18) 



with Ueff — U — 8tcos{kpa). This result is similar to 
the one obtained by Airoldi and Parola in a more general 
model . Expression (^8|) is expected to be valid when 
\Ueff\ is smaller than the effective bandwidth which is of 
the order t{Ni + Ni)/N. When this condition is satis- 
fied, we see that the model becomes equivalent to a one- 
dimensional Hubbard model with an effective interaction 



parameter Ueff- Note that the terms we have neglected 
going from ( p7| ) to (18) involve the momentum depen- 
dence of the bare interaction vertex on the external legs 
so they are irrelevant in the usual renormalization group 
analysis of weakly interacting fermions in one dimension 
prt - The most interesting consequence of ( |l8[ ) is ob- 
tained for a less than half-filled system (l/c^al < 7r/2) 
so that cos(/ci?a) > 0. In this case, the model exhibits a 
qualitative change from an effective attraction to an effec- 
tive repulsion as the bare on-site repulsion U crosses the 
value 8t cosikpa). From the previous argument, we ex- 
pect a one-dimensional Fermi hquid for U = 8tcos{kFa) 
as the off-diagonal interaction is then purely composed 
of irrelevant terms. As the density increases, the size of 
the attractive regions is reduced since the critical value of 
U decreases. Note that for dilute systems, this effective 
Hubbard model is supposed to be valid only for a narrow 
interval of values of U around 8t cos(fci?a) since the effec- 
tive single electron bandwith is small. In the dilute and 
attractive regime, it is therefore more accurate to work 
within the picture of a Bose gas of bound pairs. 



V. CONCLUSIONS 

To conclude, we have shown that a simple one- 
dimensional model in which the particle hopping is com- 
pletely assisted by the presence of another particle in 
its neighborhood generates some interesting conducting 
states at finite density. For the less than half-filled sys- 
tem, in the presence of an on-site repulsion U, we have 
established the existence of two regimes for the effective 
interaction depending on the strength of U. In the attrac- 
tive regime, the elementary building blocks of the system 
are likely to be the two-particle bound states discussed 
in section ^ Indeed, these states appear for any value 
of U and they are characterized by a very tight bind- 
ing in real space. In this regime, our model is strongly 
reminiscent of the one-dimensional attractive Hubbard 
model. An interesting difference remains the e xist ence 
of three-particle bound states detailed in section III that 
are absent in the attractive Hubbard model |Q. In the 
present work, the existence of these bound states have 
only been established in the zero density limit. It would 
be interesting to know whether they could survive in the 
presence of a finite particle density, since they could even- 
tually show up in the single electron spectral function. 

Finally, we would like to make a few comments re- 
garding the original motivation of this study, namely the 
investigation of interaction effects in tight-binding mod- 
els for which all single-particle eigenstates are localized 
in Aharonov-Bohm cages. A quasi-one-dimensional ex- 
ample illustrating this phenomenon has been analyzed 
in details in pl[ | for the two-body problem. Though 
the results obtained here have strong similarities with 
this latter model, they have also important qualitative 
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differences. A common point is the possibility to form 
extended states for the two-particle problem which are 
tightly bound in real space independently of the local 
bare repulsion strength. However, the main difference 
lies in the fact that for Aharonov-Bohm cage models, 
these extended two-particle eigenstates are, at best, de- 
generate with the localized ground states and they are 
in most cases excited states Q. In the finite density 
case, it is by no means obvious that these strongly local- 
ized systems become conducting by contrast to the one 
studied here. Nevertheless, we believe that if such con- 
ducting states appear in the presence of Aharonov-Bohm 
cages, two-particle bound states will play a major role. 
In this context, such models with strongly correlated (or 
purely assisted) hopping may provide a simple and effi- 
cient way to describe possible conducting states induced 
by interactions in otherwise localized systems. 
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